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ABSTRACT 

In the gravitational lensing of gravitational waves, the wave optics should be 
used instead of the geometrical optics when the wavelength A of the gravitational 
waves is longer than the Schwarzschild radius of the lens mass Ml- For the 
gravitational lensing of the chirp signals from the coalescence of the super massive 
black holes at the redshift zg ~ 1 relevant to LISA, the wave effects become 
important for the lens mass smaller than ~ 1O 8 M . For such cases, we compute 
how accurately we can extract the mass of the lens and the source position from 
the lensed signal. We consider two simple lens models: the point mass lens and 
the SIS (Singular Isothermal Sphere). We find that the lens mass and the source 
position can be determined within ~ 0.1% [(.S/AQ/IO 3 ] -1 for the lens mass larger 
than 1O 8 M and > 10% [(S'/iV)/10 3 ] _1 for the lens mass smaller than 1O 7 M due 
to the diffraction effect, where (S/N) is the signal to noise ratio of the unlensed 
chirp signals. For the SIS model, if the source position is outside the Einstein 
radius, only a single image exists in the geometrical optics approximation so that 
the lens parameters can not be determined. While in the wave optics cases we 
find that the lens mass can be determined even for M L < 1O 8 M . For the point 
mass lens, one can extract the lens parameters even if the source position is far 
outside the Einstein radius. As a result, the lensing cross section is an order of 
magnitude larger than that for the usual strong lensing of light. 

Subject headings: gravitational lensing - gravitational waves - binaries 



1. Introduction 

Inspirals and mergers of compact binaries are the most promising gravitational wave 
sources and will be detected by the ground based as well as the space based detectors in the 
near future (e.g. Cutler & Thorne 2002). Laser interferometers are now coming on-line or 
planned on broad frequency bands: for the high frequency band 10 — 10 4 Hz, the ground based 
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interferometers such as TAMA300, LIGO, VIRGO and GEO600 will be operated; for the low 
frequency band 10~ 4 — 10 _1 Hz, the Laser Interferometer Space Antenna 1 (LISA) will be in 
operation; for the intermediate frequency band 10~ 2 — 1 Hz, the space based interferometers 
such as DECIGO (Seto, Kawamura & Nakamura 2001) are planned. For the templates of 
the chirp signals from coalescing compact binaries, the post-Newtonian computations of the 
waveforms have been done by many authors. Using the matched filter techniques with the 
template, we can obtain the binary parameters such as the mass and the spatial position of 
the source (e.g. Cutler & Flanagan 1994). 

If the gravitational waves from coalescing binary pass near massive objects, gravitational 
lensing should occur in the same way as it does for light. The gravitational lensing of light is 
usually treated in the geometrical optics approximation, which is valid in all the observational 
situations (Schneider, Ehlers & Falco 1992; Nakamura & Deguchi 1999). However for the 
gravitational lensing of gravitational waves, the wavelength is long so that the geometrical 
optics approximation is not valid in some cases. For example, the wavelength A of the 
gravitational waves for the space interferometer is ~ 1 AU which is extremely larger than 
that of a visible light (A ~ lfi m). As shown by several authors (Ohanian 1974, Bliokh 
& Minakov 1975, Bontz & Haugan 1981, Thorne 1983, Deguchi & Watson 1986a), if the 
wavelength A is larger than the Schwarzschild radius of the lens mass M L , the diffraction 
effect is important and the magnification is small. To see the reason why the ratio M L /\ 
determines the significance of the diffraction, we consider a double slit with the slit width 
comparable to the Einstein radius £ E ~ {M L D) l l 2 where D is the distance from the screen 
to the slit (Nakamura 1998). When waves with the wavelength A pass through the slit, 
the interference pattern is produced on the screen. The width i of the central peak of the 
interference pattern is i ~ (D/^)A. Then the maximum magnification of the wave flux is 
of the order ~ (,e/^ ~ M L /X. Thus the diffraction effect is important for 



where / is the frequency of the gravitational waves. However as suggested by Ruffa (1999), 
the focused region by the gravitational lensing would have a relatively large area because of 
the diffraction, so that the lensing probability will increase. Since the gravitational waves 
from the compact binaries are coherent, the interference is also important (Mandzhos 1981, 
Ohanian 1983, Schneider & Schmid-Burgk 1985, Deguchi & Watson 1986b, Peterson & Falk 
1991). Thus we expect that the wave effects (diffraction and interference) would provide 
much information about the lens objects. 



1 See http://lisa.jpl.nasa.gov/index.html 
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In this paper, we consider the wave effects in the gravitational lensing of gravitational 
waves. We take the coalescence of the super massive black holes (SMBHs) of mass 10 4 — 
10 7 M Q as the sources. SMBH binary is one of the most promising sources for LISA and will 
be detected with very high signal to noise ratio, S/N ~ 10 3 (Bender et al. 2000). Since 
the merging SMBHs events will be detected for extremely high redshift (z > 5), the lensing 
probability is relatively high and hence some lensing events are expected. We consider the 
two simple lens models; 1) the point mass lens in which compact objects (such as black 
holes) are assumed as lens objects and 2) the SIS (Singular Isothermal Sphere) lens in which 
galaxies, star clusters and CDM (Cold Dark Matter) halos are assumed as lens. The wave 
effects become important for the lens mass 10 6 — 10 9 M Q which is determined by the LISA 
band, 10 -4 to 10 _1 Hz from Eq.(l). The frequency of the gravitational waves from the 
coalescing SMBH binary chirps so that we could see wave effects for different frequency in 
the lensed chirp signals. 

We calculate the gravitational lensed waveform using the wave optics for the two lens 
models: the point mass lens and the SIS. Then, we investigate how accurately we can extract 
the information on the lens object from the gravitational lensed signals detected by LISA 
using the Fisher-matrix formalism (e.g. Cutler & Flanagan 1994). Cutler (1998) studied the 
estimation errors for the merging SMBHs by LISA (see also Vecchio & Cutler 1998; Hughes 
2002; Moore & Hellings 2002; Hellings k Moore 2002; Seto 2002; Vecchio 2003). Following 
Cutler (1998), we calculate the estimation errors, especially for the lens mass and the source 
position. We assume the 1 yr observation before the final merging and consider the lens 
mass in the range 10 6 — 1O 9 M . Then the typical time delay between the double images is 
10 — 10 4 sec which is much smaller than 1 yr. 

This paper is organized as follows. In §2, we briefly review the wave optics in gravita- 
tional lensing for the point mass lens and the SIS model. In §3, we discuss the gravitational 
lensed waveforms detected with LISA, and mention the parameter estimation based on the 
matched filtering analysis. In §4, we numerically evaluate the signal to noise ratio and the 
parameter estimation errors. We discuss the dependence of the estimation errors on the lens 
model, the lens mass and the source position. In §5, we estimate the lensing event rate. §6 is 
devoted to summary and discussions. We assume the (Q M ,Q^) = (0.3,0.7) cosmology and 
the Hubble parameter H = 70 km/sec/Mpc, and use the units of c = G — 1. 

2. Wave Optics in Gravitational Lensing 

In this section, we briefly review the wave optics in the gravitational lensing of the grav- 
itational waves (Schneider, Ehlers & Falco 1992; Nakamura & Deguchi 1999). We consider 
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the gravitational waves propagating under the gravitational potential of the lens object. The 
metric is given by 

ds 2 = - (1 + 2U) dt 2 + (1 - 2U) dr 2 = gfjdx^dx u , (2) 

where U(r) (<C 1) is the gravitational potential of the lens object. Let us consider the linear 
perturbation h^ v in the background metric tensor g^J as 

9»» = 9$ + V- ( 3 ) 
Under the transverse traceless Lorentz gauge condition of h v = and hF 1 ^ = we have 

+ ^Z^" = °» ( 4 ) 

where ; is the covariant derivative with respect to g)iJ and R a ^ v is the background Riemann 
tensor. If the wavelength A is much smaller than the typical radius of the curvature of the 
background, we have 

W a = 0. (5) 

Following the eikonal approximation to the above equation by Baraldo, Hosoya and 
Nakamura (1999), we express the gravitational wave as 

h^ u = <f>e IJiU , (6) 

where is the the polarization tensor of the gravitational wave (e M M = 0, e^e^ = 2) and 
is a scalar. The polarization tensor is parallelly transported along the null geodesic 
(e^; a k a = 0, where k a is a wave vector) (Misner, Thorne & Wheeler 1973). Then the 
change of the polarization tensor by gravitational lensing is of the order of U (<C 1) which 
is very small in our observational situation, and hence we can regard the polarization tensor 
as a constant. Thus, we treat the scalar wave 0, instead of the gravitational wave h^ u , 
propagating through the curved space-time. The propagation equation of the scalar wave is 

^(v /Z F 5T ^ (B) ^0)=O. (7) 

For the scalar wave in the frequency domain <j>(f,r), the above equation (7) with Eq.(2) is 
rewritten as, 

(V 2 + uo 2 ) = Auo 2 U^ (8) 

where u = 2nf. The above equation (8) can be solved by using the Kirchhoff integral 
theorem (Schneider, Ehlers & Falco 1992). 

It is convenient to define the amplification factor as 

F(f) = L (/)/0(/), (9) 
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where 4> L (f) and 0(/) are the lensed and unlensed (U = in Eq.(8)) gravitational wave 
amplitudes, respectively. In Fig.l, we show the gravitational lens geometry of the source, 
the lens and the observer. Dl, D$ and D^s are the distances to the lens, the source and from 
the source to the lens, respectively. 77 is a position vector of the source in the source plane 
while £ is the impact parameter in the lens plane. We use the thin lens approximation in 
which the lens is characterized by the surface mass density £(£) and the gravitational waves 
are scattered on the thin lens plane. Then, the amplification factor F(f) at the observer is 
given by (Schneider, Ehlers & Falco 1992), 

F(f) = TfHr-- / d " x ex P [ 2 ™^( x > y)l > ( 10 ) 

where x = £/£o, and y = r}D L /^ Ds is the source position. £ is the arbitrary normalization 
constant of the length. t d is the arrival time at the observer from the source. F is normalized 
such that |F| = 1 in no lens limit (U = 0). 

Though we do not take account of the cosmological expansion in the metric, Eq.(2), 
we can apply the result without the cosmological expansion to cosmological situations since 
the wavelength of the gravitational waves is much smaller than the horizon scale. What we 
should do is 1) take the angular diameter distances and 2) replace / with f(l + z£) where 
zl is the redshift of the lens (Baraldo, Hosoya and Nakamura 1999). Then the amplification 
factor F(f) in Eq.(10) is rewritten in cosmological situations as, 

F(f) = Dsi } { \ +ZL) - [ rf 2 x exp [ 27ri/f d (x,y)] , (11) 
F> l F>ls 1 J 

where D L , D s and D LS denote the angular diameter distances. The arrival time ta at the 
observer from the source position 77 through £ is given by (Schneider, Ehlers & Falco 1992), 



1 2 

- |x - y| - V(x) + m (y) 



(12) 



The nondimensional deflection potential is determined by 

V 2 J = 2S/S cr , (13) 

where V 2 X denotes the two-dimensional Laplacian with respect to x, S is the surface mass 
density of the lens and S cr = D s / (4nD lD ls) . We choose <f> m (y) so that the minimum value 
of the arrival time is zero. We derive V ; ( x ) an d 4> m (y) m the following subsections for the 
point mass lens and the SIS model. 

In the geometrical optics limit (/ ^> t^ 1 ), the stationary points of the td(x, y) contribute 
to the integral of Eq.(ll) so that the image positions Xj are determined by the lens equation 

dt d (x,y)/dx = 0. (14) 
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This is just the Fermat's principle. The integral on the lens plane (in Eq.(ll)) is reduced to 
the sum over these images (Nakamura & Deguchi 1999) as 

F (f) = ^2 l^'| 1/2 exp t 27ri /*«W ~~ inn i\ ' ( 15 ) 

3 

where the magnification of the j-th image is fij = 1/ det (dy/dx.j), t^j = id( x j>y) and 
rij = 0,1/2,1 when Xj is a minimum, saddle, maximum point of td(x, y). In the time 
domain the wave is expressed as 

<f> L (t, r ) = Y1 ^ _ r ) exp [~ i7rn i] • ( 16 ) 

This shows that the oscillatory behavior of F(f) in high frequency / is essential to obtain 
the time delay among the images. 



2.1. Point Mass Lens 



The surface mass density is expressed as E(£) = M L 5 2 (£) where M L is the lens mass. As 
the normalization constant £ we adopt the Einstein radius given by £ = (^M L D L D LS / 'D s ) 1 ^ 2 
while the nondimensional deflection potential is ^(x) = hire. In this case, Eq.(ll) is analyt- 
ically integrated as (Peters 1974), 



F(f) = exp[— + i-(ln(-J-20 m (y)) 



% 



xr(i--» i*W-«;,l;-V], 



where w = 87iM Lz f; <f> m (y) = (x m - y) 2 /2 - lnx m with x m = (y + y/y 2 + 4)/2; M Lz = 
M L (1 + zl) is the redshifted lens mass and \Fi is the confluent hypergeometric function. 
Thus, the amplification factor F(f) includes the two lens parameters; the redshifted lens 
mass M Lz and the source position y. In the geometrical optics limit (/ M L l) from 
Eq.(15) we have 

F(f) = \fi + \ 1/2 -t\^\ 1/2 e 2mfAt «, (18) 

where the magnification of each image is fi± = 1 / 2±{y 2 + 2) / {2y \J y 2 + 4) and the time delay 
between the double images is At d = 4:M Lz [yy/y 2 + 4/2 + ln(( ^y 2 + 4 + y) /{\Jy 2 + 4 - y))]. 
The typical time delay is At d ~ AM Lz = 2 x 10 3 sec (M Lz /10 8 M Q ). 
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2.2. Singular Isothermal Sphere 

The surface density of the SIS (Singular Isothermal Sphere) is characterized by the 
velocity dispersion v as, £(£) = v 2 /(2^). As the normalization constant we adopt the 
Einstein radius £ — 4tiv 2 D l D ls /D s and the nondimensional deflection potential is if)(x.) = 
x. In this case F(f) in Eq.(ll) is expressed as 



PCX) 

F(f) = -iwe lwy2/2 / dxx J (wxy) exp 
Jo 



iw ( \<c 2 -x + (j> m (y) 



(19) 



where Jo is the Bessel function of zeroth order; <f) m (y) = y + 1/2 and w = 8n M^ z f where M^ z 
is defined as the mass inside the Einstein radius given by Ml z = 47r 2 t> 4 (l + zl)DlDls/ ' D$. 
Then, F(f) depends on the two lens parameters M Lz and y. We computed the above integral 
numerically for various parameters. In the geometrical optics limit (/ M£^), F is given 
by, 

F(f) = \fi + \ 1/2 -t\^\ 1/2 e 2mfAtd fory<l, 

= W\ 1/2 fory>l, (20) 

where pi± — ±1 + 1/y and At d = 8M Lz y. If y < 1, double images are formed. 

The wave effects in gravitational lensing of gravitational waves were discussed for the 
point mass lens (Nakamura 1998; Ruffa 1999; De Paolis et al. 2002; Zakharov & Baryshev 
2002) and a Kerr BH (Baraldo, Hosoya & Nakamura 1999). However as far as we know 
the wave effects for the SIS (Singular Isothermal Sphere) model have not been discussed 
although the SIS model can be used for more realistic lens objects such as galaxies and star 
clusters. 



2.3. The Amplification Factor 

In Fig. 2, we show the amplification factor \F(f) \ as a function of w (= %nM Lz f) for the 
fixed source position y = 0.1,0.3, 1,3 for the point mass lens (left panel) and the SIS lens 
(right panel). For w < 1, the amplification is very small due to the diffraction effect (e.g., 
Bontz & Haugan 1980). Since in this case the wave length is so long that the wave does not 
feel the existence of the lens. For w > 1, \F(f)\ asymptotically converges to the geometrical 
optics limit (Eq.(18) and (20)); 

|F(/)| 2 = \fi + \ + \fj,_\ +2|/i +/U _| 1/2 sin(27r/At,), (21) 

where //_ = for y > 1 in the SIS. The first and second terms in Eq.(21), = + 
represent the total magnification in the geometrical optics. The third term expresses the 
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interference between the double images. The oscillatory behavior (in Fig. 2) is due to this 
interference. The amplitude and the period of this oscillation are approximately equal to 
2|/x + /i_| 1 / 2 and 2irfAtd in the third term of Eq.(21), respectively. As the source position y 
increases, the total magnification (= + |/i_|) and the amplitude of the oscillation 
2|/x_|_/i_| 1//2 decrease. This is because each magnification |/i±(?/)| decreases as y increases. 
We note that even for y > 1 in SIS model (y = 3 in Fig. 2) the damped oscillatory behavior 
appears, which looks like the time delay factor of sin(27r/At rf ) although only a single image 
exists in the geometrical optics limit. 

Fig. 3 is the same as Fig. 2, but we show the phase of the amplification factor 0F{f) = 
—iln[F(f)/\F(f)\]. The behavior is similar to that of the amplitude (in Fig. 2), and the wave 
effects appear in the phase 6p as well as the amplitude |F|. For w > 1, F (f) asymptotically 
converges to the geometrical optics limit (Eq.(18) and (20)); 



$f(/) — arctan 





A*— 1 


1 /2 cos (27r/At d ) 




1/2 + 


Ia*— 


V2 sin (27r/A^)J 



(22) 



where //_ = for y > 1 in the SIS. From the above equation (22), the phase dp oscil- 
lates between — arctand/i./yU+l 1 / 2 ) and arctand/i./zi+l 1 / 2 ) with the period of 2irfAtd- As 
the source position y increases, the magnification ratio \fj,_//i + \ and the amplitude of the 
oscillation arctand/i./zi+l 1 / 2 ) decrease. 



3. Gravitational Lensed Waveform and Parameter Estimation 

3.1. Gravitational Wave Measurement with LISA 

We briefly discuss the gravitational wave measurement with LISA (see Cutler 1998; 
Bender et al. 2000). LISA consists of three spacecrafts forming an equilateral triangle and 
orbits around the Sun, trailing 20° behind the Earth. The sides of the triangle are L = 5 x 10 6 
km in length, and the plane of the triangle is inclined at 60° with respect to the ecliptic. The 
triangle rotates annually. The gravitational wave signal is reconstructed from the three data 
streams that effectively correspond to three time-varying armlength data. Two of the three 
data are linearly independent of each other. The data contain both gravitational waves 
signals to be fitted by matched filtering and noises which are assumed to be stationary, 
Gaussian and uncorrelated with each other (Cutler 1998). The gravitational wave signals 
hi j i{t) from a binary are written as 

hi A*) = ^ K>'WM*) + Fi X ,n(t)h x (t)} , (23) 
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where Fj £ (t) are the pattern functions which depend on the source's angular position 
of the binary, its orientation and detector's configuration. The quantities h+ iX (t) are the 
two polarization modes of gravitational radiation from the binary. The direction and the 
orientation of the binary and the direction of the lens are assumed to be constant during the 
observation in a fixed barycenter frame of the solar system. Further discussion and details 
about the pattern functions are shown in Cutler (1998). 



3.2. Gravitational Lensed Signal Measured by LISA 

We consider the SMBH binaries at redshift z$ as the sources. We use restricted post- 
Newtonian approximation as the in-spiral waveform (Cutler & Flanagan 1994). The coa- 
lescing time for circular orbit is typically t c = O.lyr (A^/10 6 M Q ) -5 / 3 (//10 -4 Hz)~ 8 / 3 where 
M z = (MiM 2 ) 3//5 (Mi + M 2 ) _1 / 5 (l + z s ) is the redshifted chirp mass. At the solar system 
barycenter, the unlensed waveforms h +;X (f) in the frequency domain are given by 

h + (f) = A[l + (L.n) 2 ]r^\ 

h x (f) = -2M(L.n)/- 7 /V*^, (24) 

where L (given by Ql,4>l) is the unit vector in the direction of the binary's orbital angular 
momentum and n (given by 9s, 4>s) is the unit vector toward the binary. These vectors 
are defined in a fixed barycenter frame of the solar system. The amplitude A and the 
phase ^(f) depend on six parameters; the redshifted chirp mass A4 Z and reduced mass 
/jl z = MiM 2 (l + z s )/(Mi + M 2 ); the spin-orbit coupling constant (3; a coalescence time t c 
and phase C ; the angular diameter distance to the source D S - The amplitude is 

A- fll^ML (9 * 

A -V96D S (i+z s r {2b) 

where D$(l + zs) 2 is the luminosity distance to the source, and ^f(f) is a rather complicated 
function of M Z) pL Z: (3, 4> c and t c (see Eq.(3.24) of Cutler & Flanagan 1994). 

The gravitational lensed waveforms h+ x (f) in the frequency domain are given by the 
product of the amplification factor F(f) and the unlensed waveforms h + x (/) (see section 

2); 

h L +>x (f) = F(f)h + , x (f). (26) 

where the function F(f) is given in Eq.(ll). Using Eq.(23),(24) and (26), the observed 
lensed signals h^(f) (a = I, II) with LISA are given in the stationary phase approximation 

as, 



h L M) = ^- Ds ^ 1 + ZL K f /xA a (Ht4x,y))e^M 
2 D L D LS i J 



e -i(<f>D+<P P ,a)(t+t d (x,y)) 
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xAf- 7/6 j*V\ (27) 

where <f> p , a (t) = tan- 1 [2(L • n)F a x (t)/{l + (L • n) 2 }F+(t)] and A a (f) = [(2 L • n) 2 F a x 2 (t) + 
{1 + (L ■ n) 2 } 2 F+ 2 {t)} 1 / 2 . The Doppler phase is <j> D (t) = 2nf(t)Rsm 9 S cos (0(t) - <j> s ), 
R = 1 AU and 0(t) = 27rt/T (T = 1 yr). t = t(/) is given in Eq.(3.10) of Cutler 
& Flanagan (1994). In no lens limit of V(x) = 0, the lensed signals h^(f) in Eq.(27) 
agree with the unlensed ones h a (f) in Cutler (1998). We assume the source position y is 
constant during the observation, since the characteristic scale of the interference pattern, 
~ lO 7 AU(M L2 /lO 8 M )- 1 / 2 (//mHz)' 1 [(D S D L / 'D LS ) /Gpc] 1 / 2 , is extremely larger than the 
LISA'S orbital radius (1 AU). 

Since the lensed signals h^(f) in Eq.(27) are given by double integral, we approximate 
h>a{f) in the ^ wo limiting cases; 1) geometrical optics limit (/ t^ 1 ) and 2) the time delay 
being much smaller than LISA's orbital period of (td < 1 yr). In the geometrical optics 
limit, from Eq.(15) we obtain, 

h L a (f) = ^yYl \^\ 1/2 Mt + t^)e Mft ^- im 'e- i ^ D+ ^ a ^ t+t ^ ) x Af- 7/6 e^ f \ (28) 

3 

If the time delay is much smaller than LISA's orbital period (td <C 1 yr), we expand A a ,<f)£, 
and <p Pj a around t d = as, 



~ h L a (f) = 2^A Q (f) e -Wi>+*P,«)(*) x Af' 7 '^^ 



X 



F(f) + j f {lnA a - i (d> D + p , a )} ^EMTj + o ((t,/lyr) 2 ) 



.(29) 



Since we consider the lens mass Ml z = 10 6 — 1O 9 M , the time delay is much smaller than 1 yr. 
Thus we use the above equation (29) as the lensed waveforms for the following calculations. 

In Fig. 4, the lensed signals \h^{f)\ (a = I, II) and the unlensed ones \h a (f) \ are shown. 
We show the results from one year before the final merging to the inner most stable circular 
orbit (the binary separation is r = 6(M 1 + M 2 )). We set typical parameters at the SMBH 
binary masses M li2z = 1O 6 M , the lens mass M^ z = 1O 8 M and the source position y = 1 for 
the point mass lens. The angular parameters are cos 9$ = 0.3,0s — 5.0, cos 6l = 0.8, <f>L = 
2.0, and the source redshift is z$ = 1 (the angular diameter distance is HqDs = 0.386). 
Therefore the frequency range is from 5 x 10~ 5 to 2 x 10~ 3 Hz and the time delay is 4 x 10 3 
sec. The strange behavior in the lower frequency / < 10~ 4 Hz is due to the LISA's orbital 
motion. In this frequency region, however, the difference between the lensed signal and the 
unlensed one is small due to the diffraction (see Fig. 2). On the other hand, the oscillatory 
behavior appears in the higher frequency region / > 10~ 4 Hz. This critical frequency is 



determined by the inverse of the lens mass 8%M Lz (see Fig. 2). The oscillatory amplitude and 
the period are determined by the product of the magnifications 2|/i + /i_| 1 / 2 = 2/{y^y 2 + 4) 
and the inverse of the time delay, 1/At rf (see the third term of Eq.(21)). 



3.3. 



Parameter Extraction 



We briefly mention the matched filtering analysis and the parameter estimation errors 
(Finn 1992; Cutler & Flanagan 1994). We assume that the signal h^(f) is characterized by 
some unknown parameters 7$. In the present case, there are ten source parameters (A4 Z , fi z , 
(3, 4> c , t c , D s , 9s, 05, 9~Li 4>l) and two lens parameters (M Lz ,y). In the matched filtering 
analysis the variance-covariance matrix of the parameter estimation error A'-fi is given by 
inverse of the Fisher information matrix as (A'-fiA'-fj) = (r -1 )^.. The Fisher matrix 
becomes 



where Sn(f) is the noise spectrum. The noise spectrum Sn(f) is the sum of the instrumental 
and the confusion noise, and we adopt the same noise spectrum as that in Cutler (1998). 
The signal to noise ratio (S/N) is given by 



We computed the variance-covariance matrix for a wide range of the lens parameters 
(Ml z ,i/), using the lensed waveform in Eq.(29). Since the S/N is very high for the SMBH 
merger, the Fisher matrix approach to calculate the estimation errors is valid (Cutler 1998). 
We integrate gravitational lensed waveform (in Eq.(30) and Eq.(31)) from 1 yr before the final 
merging to the cut-off frequency / cu ^ when the binary separation becomes r = 6(M 1 + M 2 ). 
We do not consider the low-frequency cut-off of LISA, which is low frequency noise wall of 
space-based instruments and is around 10~ 5 — 1CT 4 Hz (Vecchio 2003). This assumption is 
to underestimate the errors in estimation parameters. 



In this section, we present numerical results to compute the signal to noise ratio (S/N) 
and the errors in estimation parameters. We randomly distribute 100 binaries over various 
directions and orientations on celestial spheres at zs — 1 (the distance is H D S = 0.386). 
We present the mean value averaged for 100 binaries. 




(30) 




(31) 



4. 



Results 
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4.1. Lensing Effects on the Signal to Noise Ratio 

We demonstrate the gravitational lensing effect on the signal to noise ratio (S/N). In 
Fig. 5, the increasing factor of S/N by the gravitational lensing for the point mass lens is 
shown for the fixed source position y = 0.1,0.3,1,3 as a function of the lens mass M Lz . 
The vertical axis is the S/N with the gravitational lensing divided by the unlensed S/N. 
Four panels are shown for the various SMBH binary masses Mi j2z = 10 4 , 10 5 , 10 6 , 1O 7 M . 
We show the mean value averaged for 100 binaries, but the dispersion is negligibly small 
(less than 5 %). For the lens mass smaller than 10 6 M Q the magnification is very small 
irrespective of the SMBH binary masses due to the diffraction effects. In this case the 
Schwarzschild radius of the lens mass M^ z is smaller than the wavelength of gravitational 
waves A ~ 1 AU, and the waves are not magnified by lensing. This critical lens mass 
(1O 6 M ) is mainly determined by the inverse of the knee frequency of the LISA's noise 
spectrum, 1/(8tt/) ~ 8 x 10 6 M o (//mHz)" 1 (see Fig.2). But for 10 7 + 1O 7 M , the SMBH 
binary coalescences at the lower frequency (/ ~ 10~ 4 Hz), thus the critical lens mass is 
shifted for larger mass (1O 7 M ) as shown in the right bottom panel of Fig.5. This tells 
us that if the lens mass is smaller than 1O 6 M , the effect of the lens is very small. If 
the lens mass is larger than 1O 7 M , the damped oscillatory behavior appears due to the 
interference between the two images, and the S/N converges to the geometrical optics limit, 
IH 1 ^ 2 = (y 2 + 2) 1 / 2 /[y 1/ ' 2 (y 2 + 4) 1 / 4 ], which is independent of the lens mass. As y increases 
from 0.1 (solid line) to 3 (dashed line), the amplification decreases since the magnifications 
of the two images (|//±(y)|) decrease as y increases (see also Fig.2). 

Fig.6 is the same as Fig.5, but for the SIS lens model. The behavior is very similar to 
that in the point mass lens. For the lens mass larger than 1O 7 M , the S/N converges to the 
geometrical optics limit, l/i] 1 ^ 2 = (2/y) 1 ^ 2 for y < 1 and |/i| 1//2 = (1 + l/y) 1 ^ 2 for y > 1. As 
y increases from 0.1 (solid line) to 3 (dashed line), the amplification decreases (see also the 
right panel of Fig.2). 

4.2. Parameter Estimation for the Lens Objects 

In this section, we show the parameter estimation for the lens objects. We show the 
results for the SMBH binary with masses 10 6 + 1O 6 M , because we found S/N is higher than 
the other binary masses (M lj2z = 10 4 , 10 5 and 1O 7 M ). We distribute the 100 binaries over 
the various directions and the orientations at zs — 1, and the mean value of the S/N without 
lensing is 2600 in these 100 binaries. We show the mean value of errors averaged for 100 
binaries, for M Lz < 1O 7 M the dispersion is relatively large (< 40%), but for M Lz > 1O 7 M 
the results converge to that in the geometrical optics limit and the dispersion is negligibly 
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small. 

In Fig. 7, the estimation errors for the redshifted lens mass AM Lz (left panel) and the 
source position Ay (right panel) are shown as a function of M Lz with y = 0.1,0.3, 1,3 for 
the point mass lens. We use the units of S/N = 10 3 , and the results (AM Lz , Ay) scale 
as (S/N) -1 . For M Lz < 10 7 M© the estimation errors are relatively large > 10%, since the 
effect of lensing on the signals is very small due to the diffraction. For M Lz > 1O 8 M the 
geometrical optics approximation is valid, and the errors converge to a constant in Fig. 7. 
The redshifted lens mass and the source position can be determined up to the accuracy of 
~ 0.1%, as shown in Fig. 7. The errors in the geometrical optics limit are well fitted by (see 
Appendix A), 



AM Lz 1 , y/y(y 2 + 2)(y 2 + 4) 5 /* 

M Lz ~ S/N* 2r 

V ~ S/N X 2^y ' [d2} 

where S/N is in the unlensed case, and r = At d /AM Lz = y^/y 2 + A/2 + ln((vV + 4 + 
y)l (\fy 2 + 4 — y)). Thus, one could determine the lens parameters, the redshifted lens mass 
and the source position, up to the accuracy of ~ (S/N) -1 . The above equations (32) are 
valid if the time delay Atd is much smaller than the LISA's orbital period 1 yr. If the time 
delay Atd becomes comparable to 1 yr, the LISA's orbital motion affects the results. 

Fig.8 is the same as Fig. 7, but as a function of y. For y > 1, the errors are convergent 
to the geometrical optics limit of Eq.(32) irrespective of the lens mass. As y increases, the 
time delay t& increases, and the geometrical optics limit (ft^ ^> 1) is valid. We note that 
even for y > 10 one can extract the lens information. In the case of light, the observable 
is the lensed flux which is proportional to the magnifications, oc \fJ>±\, but for gravitational 
waves the observable is the lensed amplitude which is proportional to the square root of the 
magnifications oc |/i±| 1//2 . For example, let us consider the case where the flux ratio of a 
brighter image to a fainter one is 100 : 1. Then the amplitude ratio is 10 : 1 so that the 
fainter image can be observed even if the source position is far from the Einstein radius 
in the case of gravitational waves. Denoting the largest source position for which one can 
extract the lens parameters as y cr , we approximate the errors in Eq.(32) for the large y limit; 
A7/7 ~ (S/N^y 2 , where 7 = M Lz ,y. Then we obtain 

Thus the lensing cross section (oc y 2 cr ) increases an order of magnitude larger than that for 
the usual strong lensing of light (y cr = 1) (e.g. Turner, Ostriker & Gott 1984). 
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In Fig. 9, the estimation errors for the SIS model are shown. For M Lz < 1O 7 M the 
behavior is similar for the point mass lens. But for M^ z > 1O 8 M the behavior strongly 
depends on y. In the geometrical optics approximation, the errors are given by (Appendix 
A), 

^ = ^ f 0ry <l. (34) 

M Lz y S/N y y y ~ K J 

and the lens parameters are not determined for y > 1. We note that even for y = 3 the 
lens parameters can be extracted for M Lz ~ 10 6 - 1O 8 M due to the wave effects. For 
y = 0.1 and 0.3, the asymptotic behavior of errors are somewhat smaller than the results 
in Eq.(34), because the order of 1/f term in F(f) (which is neglected in the geometrical 
optics approximation / — > oo) affects the results. For y — 1, the errors decrease as lens mass 
increases as shown in Fig. 9, because the errors converge to the results in the geometrical 
optics limit of Eq.(34) which vanish at y — 1. As a result, if y > 1, the errors asymptotically 
increase with the increase of the lens mass, but if y < 1, they asymptotically converge to 
constants. 

Fig. 10 is the same as Fig.9, but as a function of y. We note that even for larger y > 1 
we can extract the lens information. Thus the lensing probability (oc y 2 ) to determine the 
lens parameters increases as compared with the results in geometrical optics limit for the 
lens objects in the mass range 10 6 — 1O 8 M . 



4.3. Lensing Effects on the Estimation Errors of the Binary Parameters 

We discuss the gravitational lensing effects on the estimation errors of the SMBH binary 
parameters. We study five binary parameters; the redshifted chirp mass A4 CZ , the reduced 
mass fj, z , the distance to the source D s and the angular resolution (9$, 4>s)- We find that the 
estimation errors of these parameters decrease because S/N increases by lensing (see Fig. 5 
and 6). The error A7 is roughly proportional to the inverse of the S/N as A7 oc (S/N)^ 1 
(see Eq.(30) and (31)). 



4.4. Results for Various SMBH Masses and Redshifts 

So far we presented the results for equal mass SMBH binaries with redshift zs — 1. 
In this section, we comment the results for the case of various (unequal) SMBH masses 
10 4 - 1O 7 M and redshifts z s = 1 - 10. 

The critical lens mass in which the wave effects become important (10 6 — 1O 8 M ) is 
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mainly determined by the inverse of the knee frequency of the LISA's noise spectrum, ~ 
8 x 1O 6 M (//mHz)" 1 , independent of the binary mass (see section 4.1). But for the massive 
total mass binary (M l2 + M 2z ) > 1O 7 M , the binary coalescences at the lower frequency 
(~ 10~ 4 Hz), thus the critical lens mass is shifted for larger mass (10 7 — 1O 9 M ). For the larger 
lens mass Ml z <) 10 8 M Q , the results (the S/N increase and the estimation errors) converge 
to that in the geometrical optics limit irrespective of the binary mass. The estimation errors 
in Fig. 7-10 are the case of 10 6 — 10 6 M Q binary at redshift z s = I and are normalized to 
S/N = 10 3 and simply scale as (S/N)^ 1 . In order to translate the results in the various 
unequal SMBH binaries, we present the S/N for binary masses Mi j2z = 10 4 — 10 7 M Q with 
redshifts z s = 1, 3, 5, 10 in Table. 1. We assume 1 yr observation of in-spiral phase before final 
merging. The results are the mean value of 100 binaries which are randomly distributed at 
each redshift, and the dispersion is relatively large ~ 50%. From Table. 1, one could translate 
the results in Fig. 7-10 into errors in real situations. 

We also comment the results for the case of only hi data available, while we used the 
combination of hj and hu data (see section 3.1). In this case, the S/N increase in Fig. 5 and 
6 are not changed, but the estimation errors are slightly larger (~ 30%) than that in Fig. 7-9 
for Ml z < 1O 7 M if the errors are normalized to S/N = 10 3 . We note that the S/N is \pl 
times smaller than that in the case of the two data available in Table. 1. 



5. Lensing Event Rate 

We discuss the event rate of merging SMBHs and estimate the lensing probability and 
the lensing event rate. The expected rate of merging SMBHs detected by LISA is in the 
range 0.1 - 10 2 events per year (Haehnelt 1994,1998). Recently, Wyithe & Loeb (2002) 
suggested that some hundreds detectable events per year could be expected, considering the 
merger rate at exceedingly high redshift (z > 5 — 10). Thus we take ~ 300 events per year 
as the merging event rate. 

We consider the lens objects distributed over the universe and calculate the lensing 
probability for each lens model. For the point mass lens, we take the compact objects 
(10 6 — 1O 9 M ) such as black holes as lens. Denoting the mass density parameter of compact 
objects as f2 C o, the lensing probability for a source at redshift zs is (Schneider, Ehlers & 
Falco 1992), 

where H(z) is the Hubble parameter at redshift z. The cosmo logical abundance of the 
compact objects in the mass range 10 6 — 1O 9 M is limited by f2 C o < 0.01 by the search for 
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multiple images in radio sources (Wilkinson et al. 2001; see also Nemiroff et al. 2001). In 
Table 2, we show the upper limit on the lensing probability for the point mass lens. Since we 
set y cr = 10 (Eq.(33)), the lensing probability is one hundred times larger than that normally 
assumed for the strong lensing of light (y cr = 1). As shown in Table 2, the upper limit of 
the lensing probability is very high (almost 1) and is typically ~ (f2co/10~ 2 ). The lensing 
event rate is the product of the merging rate (~ 300 per year) and the lensing probability, 
so that the lensing events will be 1 event per year if fi co = 10~ 4 . 

For the SIS model we take CDM halos (10 6 — 10 9 M Q ) as the lens objects (e.g. Narayan 
& White 1988). The lensing probability is 

p , x 2 f' 3 , (1 + *L? H D LS (z L ,z s ) H D L (z L ) /- 1q9m ° 

P(z s ) = ny cr dz L / dM L vN v (v,z L ), (36) 

Jo n{Z L )/H HqD S {Z S ) JlO6M 

where N v is the comoving number density of the lens and is assumed to be given by the 
Press-Schechter velocity function (Press & Schechter 1974) with a s = 1. In Table 2, we show 
the lensing probability for the SIS model. We set y cr = 3 (see Fig. 9) and hence the lensing 
probability is almost ten times larger than that for light (y cr = 1). As shown in Table 2, the 
lensing probability is typically ~ 10~ 4 — 10~ 3 . The merger rate is ~ 300 events per year at 
high redshift (z > 5), then the lensing events would be 1 event per year. 

We note that the results in Table. 2 are for the case of the S/N = 10 3 , and are somewhat 
overestimated for the binaries of S/N < 10 3 in Table. 1. For example, the lensing probability 
is proportional to (S/N) from Eq.(33) for the point mass lens, and it is appropriate to use 
y cr = 1 for S/N < 10 3 in the SIS. In the case of the high event rate (~ 300events/year), 
many fainter signals (S/N <C 10 3 ) are expected and we note that the errors in Fig. 7-10 are 
worse for these binaries. 

Next, we discuss how we can identify the lensing signal. If the lensing event occurs, the 
amplitude and the arrival time of the gravitational waves are changed by lensing. But the 
other features (such as binary mass) are not changed. Thus, if the two signals have the same 
binary parameters (such as chirp mass) except for the amplitude and the arrival time, that 
would be a signature of gravitational lensing in the geometrical optics limit. More generally, 
oscillatory behavior in the waveform \h L (f) \ is a signature of gravitational lensing (see Fig. 4). 
However it will be difficult to identify the source and the lens objects in the sky, since the 
angular resolution of the LISA is ~ 1 deg (see Cutler 1998). Furthermore, the gravitational 
wave amplitude is changed by the lensing magnification and hence one must assume the lens 
model in order to determine the distance to the source. (Effect of lensing on measuring the 
distance is recently discussed in Holz & Hughes (2002).) As one determines the distance to 
the source D s (zs), the redshift z s (D s ) could be determined if the cosmological parameters 
are well known (see Hughes 2002). 
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6. Summary 

We have discussed the gravitational lensing of gravitational waves from chirping binaries, 
taking account of the wave effects in gravitational lensing. The SMBH binary is taken as 
the source detected by LISA, and the two simple lens models are considered: the point 
mass lens and the SIS model. We calculate the lensing effects on the signal to noise ratio 
(S/N) and how accurately the information of the lens object, its mass, can be extracted 
from the lensed signal. As expected, for the lens mass smaller than 10 8 M Q , the wave effects 
are very important to calculate the S/N and the errors in the estimation parameters. It is 
found for the lens mass smaller than 10 6 M Q the signals are not magnified by lensing due 
to the diffraction effect. For the lens mass larger than 1O 8 M the lens parameters can be 
determined within (very roughly) ~ 0.1% [(S / N) /10 3 ]" 1 . We note that the lensing cross 
section to determine the lens parameters is order of magnitude larger than that for light. 

In this paper, we calculate the case for LISA. But similar analysis can be done for other 
detectors. For the ground-based interferometers (TAMA300, LIGO, VIRGO, GEO600), 
neutron star binaries are taken as the sources and the lens mass for which the wave effects 
become important is 10 — 10 4 M Q . Similarly for the space-based interferometers such as 
DECIGO (Seto, Kawamura & Nakamura 2001), the important lens mass becomes 10 5 — 
1O 7 M . Since mergers of neutron star binaries will be detected at least several per year for 
LIGO II (Phinney 1991; Kalogera et al. 2001) and ~ 10 5 per year for DECIGO, the lensing 
events would also be expected for other detectors. 

We would like to thank Naoki Seto and Takeshi Chiba for useful comments and dis- 
cussions. This work was supported in part by Grant-in-Aid for Scientific Research of the 
Japanese Ministry of Education, Culture, Sports, Science and Technology, No. 140472 12 
(TN), and No. 14204024 (TN). 

A. Estimation errors in the geometrical optics limit 

To evaluate the estimation errors AM Lz , Ay in the geometrical optics limit, we consider 
the simple waveform; 

^(/) = (k| 1/2 -^K| 1/2 e MAti )xM/), (Al) 

where h oc A is the unlenesd signal and Atd oc Ml z , with three parameters ji = (lnM^y, ln„4). 
Then, the Fisher matrix Tij(i,j = 1,2,3) in Eq.(30) can be analytically obtained as, 

r n = (27rAt d ) 2 M (fh\fh), 
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The S/N is the signal to noise ratio for the unlensed signal h. The estimation errors can 
be analytically obtained by the inverse of the Fisher matrix, AM Lz /M Lz = [(r^n] 1 / 2 and 
Ay/y = [(r _1 ) 2 2] 1//2 /y- Using the geometrical optics approximation, fAt d 3> 1, we obtain 
the errors with Eq.(A2) and (A3) as 
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The above equations (A4) are used for the general lens model when the double images form. 
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Fig. 1. — Gravitational lens geometry for the source, the lens and the observer. Dl,D$ 
and Dls are the distances between them. 77 is a displacement of the source and £ is an 
impact parameter. We use the thin lens approximation in which the gravitational waves are 
scattered in the thin lens plane. 
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Fig. 2. — The amplification factor \F(f) \ as a function of w (= SnMi z f) with the fixed source 
position y = 0.1,0.3, 1,3 for the point mass lens (left panel) and the SIS (right panel). For 
w < 1, the amplification is very small due to the diffraction effect. For w > 1, the oscillatory 
behavior appears due to the interference between the double images. We note for the SIS 
that even if y > 1 (a single image is formed in the geometrical optics limit) the damped 
oscillatory behavior appears. 
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Fig. 3. — Same as Fig. 2, but the phase of the amplification factor ^f(/) = — i l n [-^(/)/|-^(/)|] 
as a function of w (= 8nM Lz f). 
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Fig. 4. — The lensed signals \h^(f)\ (a = I, II) (solid line) and unlensed ones \h a (f)\ 
(dashed line) measured by LISA. The signals are shown from lyr before coalescence to 
ISCO of r = 6(Mi + M 2 ). The redshifted masses of the SMBH binary is M h2z = 1O 6 M , the 
redshifted lens mass is M Lz = 1O 8 M and the source position y — 1. The angular parameters 
are cos 8s = 0.3, <ps = 5.0, cos$l = 0.8, <pL = 2.0, and the source redshift is z$ = 1 (distance 
is H Ds = 0.386). The strange behavior for / < 10 -4 Hz is due to the LISA orbital motion, 
and the difference between the two signals is small due to the diffraction. On the other hand, 
the oscillatory behavior appears for / > 10~ 4 Hz which is determined by the inverse of the 
lens mass 8irM Lz (see Fig. 2). This oscillation is due to the interference between the double 
images. 
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Fig. 5. — The increasing factor in S/N due to the gravitational lensing by the point mass lens 
for the various SMBH binary masses M h2z = 10 4 , 10 5 , 10 6 , 1O 7 M . The horizontal axis is the 
redshifted lens mass; the vertical axis is the lensed S/N divided by the unlensed S/N. The 
source position is fixed at y — 0.1, 0.3, 1, 3. For Ml z < 1O 6 M , the magnification is very small 
due to the diffraction effect irrespective of the SMBH binary masses. For Ml z ^ 1O 7 M , 
the damped oscillatory patterns appear due to the interference between the two images, and 
this behavior converge in the geometrical optics limit, l/i] 1 ^ 2 = (y 2 + 2) 1 / 2 /[y l ^ 2 {y 2 + 4) 1 / 4 ]. 
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Fig. 6. — Same as Fig.5, but for the SIS lens model. For Ml z > 1O 7 M , the results converge 
in the geometrical optics limit, = (2/y) 1 / 2 for y < 1 and l^ 1 / 2 = (1 + 1/t/) 1 / 2 for y > 1. 
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Fig. 7. — The estimation errors for the redshifted lens mass AM^ 2 (left panel) and the source 
position Ay (right panel) for the point mass lens. The results are presented for the SMBH 
binary of masses 10 6 + 1O 6 M at zs = 1. The errors are normalized by S/N = 10 3 and 
simply scale as (S/N)^ 1 . For M^ z < 10 7 M Q the errors are relatively large, since the effect 
of lensing is very small due to the diffraction. For M Lz > 1O 8 M the geometrical optics 
approximation is valid, and errors converge to constants. 




Fig. 8. — Same as Fig. 7, but as a function of y. We note that even for y > 10 we can extract 
the lens information. Then the lensing cross section (oc y 2 ) increases an order of magnitude 
larger than that for the usual strong lensing of light (y — 1). 
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Fig. 9. — The estimation errors for the redshifted lens mass AM^ 2 (left) and the source 
position Ay (right) for the SIS model. The results are presented for the SMBH binary of 
masses 10 6 + 1O 6 M at z$ = 1. The errors are normalized by S/N = 10 3 and simply scale 
as (S/N)^ 1 . Even for y = 3 (a single image is formed in the geometrical optics limit), the 
lens parameters can be extracted at M Lz ~ 10 6 — 1O 8 M due to the wave effects. 




Fig. 10. — Same as Fig.9, but as a function of y. We note that even for y > 1 we can 
determine the lens parameters. Then the lensing cross section (oc y 2 ) becomes larger than 
that in the geometrical optics approximation (y — 1). 
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Table 1: The signal to noise ratio (S/N) for the various binary masses 10 4 — 10 7 M Q with 
redshift Zs = 1, 3, 5, 10. We assume 1 yr observation of in-spiral phase before final merging. 



Lens Model 


z s = l 


^5 = 3 


z s = 5 


as = 10 


Point mass lens 


< 0.21 


< 1.1 


< 2.0 


< 3.9 


SIS 


7.2 x 10~ 5 


8.1 x 10" 4 


2.0 x 10" 3 


4.7 x 10~ 3 



Table 2: The lensing probability by the lens mass in the range 10 6 — 1O 9 M with the source 
redshift zs = 1, 3, 5, 10. For the point mass lens, we give the upper limit which is determined 
by the observational constraint on the abundance of the compact objects. When the lensing 
probability is more than one, the lensing occurs some times. For the SIS, CDM halos are 
assumed to be lenses. The presented values are for the case of S/N — 10 3 , and hence the 
results are somewhat overestimated for the binaries of S/N < 10 3 in Table. 1. If the expected 
rate of merging SMBHs is ~ 300 per year (Wyithe & Loeb 2002), then the lensing events 
will be detected 1 event per year. 



